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Global and local level dynamics in chaotic microwave 
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Abstract. The spectra of Sinai microwave billiards and rectangular billiards with statisti- 
cally distributed circular scatterers have been taken as a function of the position of one wall, 
and of one of the scatterers, respectively. Whereas in the first case the velocity distribution 
and correlations obey the universal behaviour predicted by Simons and Altshuler, in the sec- 
ond case a completely different behaviour is observed. This is due to the fact that a shift of 
one wall changes the wave function globally whereas the displacement of one scatterer only 
leads to a local perturbance. 
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1 Introduction 

The eigenvalue dynamics of the spectra of disordered systems in dependence of an 
external parameter has attracted considerable interest in recent years. It has been 
motivated by the fact that the conductance through a system is closely linked to the 
sensitivity of its eigenvalues on a perturbation ||l|, ||. In this respect the work of 
Akkermans and Montambaux has to be mentioned in particular |^ , who showed that 
the conductance can be expressed in terms of the quadratically averaged eigenvalue 
velocities, 
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Here the brackets denote a local average, and X is the external parameter. Later it 
was shown by Simons and Altshuler using supersymmetry techniques that the eigen- 
value velocities v should be Gaussian distributed Q 
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Moreover, if the eigenvalues En and the level dynamics parameter X are rescaled 
according to 
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Fig. 1 Sketch of the quarter Sinai biUiard used for the global level dynamics (left) and of 
the rectangular billiard with randomly distributed scatterers (right). Both figures are drawn 
in scale, including size and positions of the scatterers. For the dimensions see the text. 



where A is the mean level spacing, the velocity autocorrelation function 

, , _ / denjxQ + x) denjxo) \ _ / denix^V^ 
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should obey a universal behaviour |Q, ||] . A recent review on these questions can be 
found in Q. 

There are a number of attempts to study the autocorrelation function (^), both 
theoretically and experimentally. The systems under consideration ranged from the 
hydrogen atom in a strong magnetic field |^ via conformally deformed |^ and ray- 
splitting billiards |^ to the acoustic spectra of vibrating quartz blocks Q . In all cases 
the overall behaviour predicted by Simons and Altshuler was more or less observed 
but in nearly all cases there had been significant deviations as well. 

This was our motivation to study velocity distributions and correlations in various 
mircowave billiards by applying different types of parameter variations. Most to our 
surprise we found that the velocity distributions show a clear dependence on the 
type of parameter variation. Only for the type classified as 'global' below, we found 
Simons' and Altshuler's universal behaviour, whereas for a second type denoted 'local' 
a completely different behaviour was found. A preliminary version of this paper has 
been published already in p| . 



2 Experiment 

All results presented below have been obtained by taking the eigenfrequency spectra of 
billiard-shaped microwave cavities, shortly termed microwave billiards. The technique 
is described in jl^ . Here it may be sufficient to note that in quasi- two-dimensional 
resonators, i.e. cavities with top and bottom plate parallel to each other, quantum 
mechanical and electromagnetic spectra are completely equivalent, as long as the max- 
imum frequency z^max = c/2h is not surpassed, where h is the height of the resonator. 
In the present experiment the height was 8 mm yielding a maximum frequency of 18.74 
GHz. 
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Fig. 2 Part of the eigenvalue level dynamics spectra for the quarter Sinai billiard as a 
function of the length (left) and the rectangular billiard with randomly distributed scatterers 
as a function of the position of one scatterers (right). All spectra have been unfolded to a 
mean density of one. 



In the very beginning it was not evident at all which billiard parameter should be 
varied. According to Simons and Altshuler the parameter X 'could denote the strength 
of some field, like the Aharonov-Bohm flux through a ring, or a magnetic field, or the 
position of some impurity' (cited literally from Ref. ||]). In microwave billiards the 
magnetic field is not available as a parameter (though the use of billiards with ferrite- 
coated walls may offer an alternative (l^). On the other hand it is straightforward 
to study the spectra of microwave billiards as a function of some length The 
magnetic field and the length variation have in common that already a small change 
in the parameter will modify the wavefunction everywhere. For this reason we call 
this type of level dynamics 'global'. The variation of the position of one impurity, on 
the other hand, can be easily performed in microwave billiards, with the additional 
advantage that, in contrast to real mesoscopic systems, we can control the disorder 
in every moment perfectly. The displacement of one impurity changes the value of 
the wave function only in its neighbourhood. Consequently, we call this type of level 
dynamics 'local'. 

Both types of level dynamics were realized in the experiments. One of the systems 
studied was a quarter Sinai billiard with width 6=200 mm, radius r=70 mm of the 
quarter circle, and a length a which was varied between 460 and 480 mm in steps of 
0.2 mm. About 120 eigenvalues entered into the data analysis in the frequency range 
14.5 to 15.5 GHz. The second system was a rectangular billiard with side lengths 
a=340 mm, &=240 mm, containing 19 circular disks with a diameter of 4.6 mm whose 
positions were determined by means of a random number generator. The position of 
one further disk was varied in one direction in steps of 1 mm. Two different movable 
disks were used with diameters of 4.6 and 20 mm, respectively (see Fig. |l|). 

A part of the level dynamic spectrum of the Sinai billiard is shown in Fig. |^ 
(left). As a function of length the resonances experience four to five avoided crossing 
in the studied length range. This is about the maximum range usable for the velocity 
correlation measurements for the following reason. Upon variation of the length it hap- 
pens that a node line of an eigenfunction passes the position of the coupling antenna. 
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Fig. 3 Velocity distribution for the quarter Sinai billiard on a linear (left) and a logarithmic 
(right) scale. The dotted line corresponds to a Gaussian, the dashed line to a modified Bessel 
function (see Eq. (^)). 



leading to a temporary loss of the eigenfrequency. In the studied length range each 
resonance disappeared at most once, but with increasing length range the tieing up of 
the open ends became more and more unreliable. The situation changes qualitatively 
of the level dynamics of the rectangular billiard with randomly distributed scatterers 
shown in Fig. |^ (right). Already a comparison by eyes shows that there is much more 
regularity in the latter case. Now it is no problem to follow each eigenvalue over about 
10 avoided crossings. 



3 Global level dynamics 

We start with a discussion of the Sinai billiard results. Fig. ^ shows the found velocity 
distribution both on a linear and a logarithmic scale. The found distribution is close to 
a Gaussian in accordance with the universal prediction. The dashed line is a modified 
Bessel function which will become important in the following. It is well-known that 
non-generic features such as bouncing balls give rise to deviations from a Gaussian 
distribution in particular in the region of the wings |l^ . In the measurement we 
therefore carefully avoided the regions disturbed by the dominating bouncing ball. 
Hence such a disturbance does not show up in the present case. 

Fig. ^ shows the corresponding velocity autocorrelation function. To obtain the 
result, each eigenvalue was followed over the available length range including four 
to five avoided crossings. The mean squared velocity needed for the scaling was 
calculated for each eigenvalue independently. Subsequently the results of all eigenval- 
ues were superimposed. Such a procedure was mandatory, since the individual mean 
squared velocities showed large variations. For this reason a scaling with a globally 
calculated mean squared velocity did not give satisfactory results. 

The dotted line in Fig. ^ corresponds to Simons' and Altshuler's universal velocity 
correlation function. Since an analytical expression for this quantity is not available, 
it was obtained from a random matrix simulation [ p^ . The overall agreement between 
experiment and theory is good, apart from some small but significant deviations close 
to the minimum. It was already mentioned in the beginning that in most studied cases 
such an agreement was not found. 
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Fig. 4 Velocity autocorrelation function for the quarter Sinai billiard as a function of the 
position of one length rescaled according to Eq. (^. The dotted line corresponds to Simons' 
and Altshuler's universal function. 

4 Local level dynamics 

Whether a movable scatterer perturbs the system globally or locally, depends on its 
size, or, more precisely, on the parameter 6 — kD, where D is the diameter of the 
scatterer and k the wavenumber. For small S values ((5 ^ 1) a local level dynamics 
is expected, whereas with increasing S values a gradual transition towards a global 
behaviour should be observed. Fig. || shows the velocity distribution in the rectangular 
billiard with randomly distributed scatterers for three different 6 ranges, again in a 
linear and a logarithmic plot. In Fig. |^(a) the diameter of the movable disk was 4.6 
mm, and the eigenvalues were taken in the range 3.4 to 6 GHz. In Figs, ^(b) and (c) 
the diameter of the scatterer was 20 mm, and the eigenvalues were in the range 3.4 to 
5.8 GHz and 12.5 to 14.3 GHz, respectively. 

For small S values the found velocity distribution is well described by a modified 
Bessel function 



The logarithmic plot shows that this behaviour is found over two orders of magni- 
tude. Note that now asymptotically a single exponential behaviour is found in contrast 
to the Gaussian behaviour observed for the global velocities. 

With increasing 6 values a transition to Gaussian behaviour is observed as expected. 
We completed the sequence with a level dynamics in a Sinai billiard, where the central 
circle was moved, corresponding to S values between 30 and 40. Here the found 
velocities were again Gaussian distributed (not shown). 

A corresponding deviation from the universal behaviour is found for the velocity 
autocorrelation functions. The left column of Fig. ^ shows the results for the same 6 
ranges as in Fig. ^. There is no similarity at all with the universal behaviour. Only 
for the largest S range the found correlation function seems to approach the universal 
function. 

All these findings are in striking contrast to all what has been expected by the 
experts. Nevertheless, the discrepancy from the universal behaviour can be quite easily 
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Fig. 5 Velocity distributions in a rectangular billiard with randomly distributed scatterers 
as a function of the position of one of the scatterers on a linear (left column) and a logarithmic 
(right column) scale. The histograms were taken in three different ranges of S = kD, where D 
is the diameter of the movable scatterer: 0.35 < 5 < 0.65 (a), 1.4 < 6 < 2.6 (b), 5.1 < S < 5.9 
(c). The dotted line corresponds to a Gaussian, the dashed line to a modified Bessel function. 



resolved. It is a well-known fact that the insertion of a metallic object into a quasi- 
two-dimensional microwave cavity shifts the eigenfrequencies to smaller frequencies. 
If the dimensions of the perturber are small compared to the wavelength, this shift is 
proportional to the square of the electric field. This technique has been used routinely 
for many years to map field distributions in microwave cavities JT7[ and has recently 
been applied to the study of wave function in chaotic billiards as well ^ pO[ . 
Applied to the present problem this means that in the limit S ^ the eigenvalue 
velocities are given by 

|| = «V|V^P (6) 

where X is the position of the scatterer, and V is the gradient taken in the direction 
of the displacement, a is a constant factor depending exclusively on the geometry of 
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Fig. 6 Velocity autocorrelation functions in a rectangular billiard with randomly distributed 
scatterers for the same ranges of S as in Fig. |^. On the left column the position of the movable 
scatterer was scaled according to Eq. (^). The dotted line is Simons' and Altshuler's universal 
function. On the right column the same data are plotted, but with a scaling of the position 
according to Eq. (fUj). The dashed lines have been calculated from Eq. (p|). 



the perturber. It follows for the velocity distribution 

P{v) = {d{v - 2ai/)VV')) (7) 

where the brackets denote an average of eigenvalues. This average can be calculated, 
if we assume that the wavefunction ip{r) can be described by a random superposition 
of plane waves |^ 

V'M-^a^.e*^^"^^ (8) 

where the modulus of the wavenumbers of all waves contributing is the same, |A:„| = 
k, but where amplitudes a„ and directions kn/k are random. This approach has 
been proved extremely fruitful in the description of eigenfunctions in chaotic systems 
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By means of the central limit theorem it is now an easy matter to show that under 
these assmptions P{ip) and PlV-tli) are uncorrelated and both Gaussian distributed. 
With these ingredients the average (|^) can be calculated, resulting in the modified 
Bessel function (||) for the velocity distribution. With help of the ansatz (|) the 
velocity autocorrelation can be calculated as well, using standard techniques |24|. The 
result is 

c{x) = - [Jlix)]" = Jl{x) - 2Jl{x) - Mx)J2{,x) (9) 

where 

X = kX. (10) 

Note that in contrast to Simons' and Altshuler's approach where the level dynamics 
parameter X has been scaled by means of the square root of the quadratically averaged 
velocities, now we have got a straightforward scaling in terms of the wavenumber k. 

On the right column of Fig. ^ the same data shown on the left column are plotted 
again, but now with an abscissa scaled according to Eq. (p^). The solid line has 
been calculated from Eq. (^). It describes the experimental results quite well. In 
particular the wavelengths of the oscillations are reproduced correctly even for the 
larger 6 ranges. Only an increasing damping of the oscillations with increasing S is 
observed. This oscillatory behaviour in the velocity autocorrelation function, which 
is already evident from a visual inspection of Fig. || (right), is in sharp contrast to 
the global case, where the autocorrelation function shows only one minimum, while 
all subsequent oscillations found for the individual levels are completely wiped out by 
averaging process. 



5 Conclusion 

This work has shown that in contrast to our previous understanding there is not one 
but there are two regimes of universality. Only for the global case where a small 
variation of the level dynamics parameter changes the wave functions everywhere, the 
universal velocity distributions and correlations predicted by Simons and Altshuler are 
really observed. In this range the scaling of the level dynamics parameter in terms of 
the square root of the quadratically averaged velocities is the appropriate one. 

In the local limit, on the other hand, realized in the present work by moving a small 
pertuber through a disordered system, another universal regime is observed. Here the 
scaling of the level dynamics parameter in terms of the wave number is the correct 
one, and another class of velocity distributions and correlations is found, which has 
been discussed in this paper. 

The issue of this paper has been discussed with numerous colleagues from the mesoscopic 
community. Discussions with Y. Fyodorov and T. Guhr have been in particular helpful. E. 
Mucciolo provided us with his unpublished calculation of the universal velocity autocorrela- 
tion function. The work has been supported by the Deutsche Forschungsgemeinschaft via the 
SFB 185 "Nichthneare Dynamik" as well as by an individual grant. 
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